Molecular orbital excitations in cuprates 
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We report resonant inelastic x-ray scattering studies of electronic excitations in a wide variety 
of cuprate compounds. Specifically, we focus on the charge-transfer type excitation of an electron 
from a bonding molecular orbital to an antibonding molecular orbital in a copper oxygen plaquette. 
Both the excitation energy and the amount of dispersion are found to increase significantly as the 
copper oxygen bond-length is reduced. We also find that the estimated bond-length dependence of 
the hopping integral tpd is much stronger than that expected from tight-binding theory. 
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I. INTRODUCTION 

One of the most important characteristics of the elec- 
tronic structure of the cuprates is the strong hybridiza- 
tion between the Cu ?id,j.2_y2 level and the O 2p^ level, 
where Pr, denotes the Px or py orbitals pointing towards 
the Cu ions. Because of this hybridization, the Cu-0 
bond has a strong covalent character and a large energy 
splitting exists between the bonding (a) and antibond- 
ing (cr*) molecular orbitals. In the ionic limit without 
hybridization, this energy splitting corresponds to the 
energy difference between the atomic ci^2_y2 and Pa or- 
bitals, which is Aq ~ 3.5 eV.^ As the p — d hybridization 
becomes larger, the energy splitting (Ao-o-* ) between the 
two molecular orbitals increases, reflecting the increas- 
ingly covalent nature of the Cu-0 bonding. Thus, Ao-o-* 
is a direct measure of the Cu-0 hybridization, and could 
serve as an independent route to determine the value 
of the hopping matrix element, tpd, since it is believed 
that Ao-o-* is directly related to tpd^ Although there 
has been no systematic theoretical study of Ao-o-* , the 
values obtained from first-principles calculations range 
widely from ~ 4 eV to ~ 9 eVA^i^ For example, Martin 
and Hay^ have carried out an ah initio quantum chem- 
istry calculation of a cluster of copper oxygen octahedra 
(CuO^°-) in Lai. 85Sro.i5Cu04, and obtained l^aa" ^ 9 
eV, while in a recent density-functional calculation of a 
similar cluster, Hiisser and coworkers reported a value of 
Actct* ^ 5.8 eVi^ Experimentally, determining this quan- 



tity has been very difficult, since in this energy range 
transitions involving the La 4/ bands dominate the spec- 
tral features of optical spectroscopy. 

In this paper, we report a systematic experimen- 
tal study of Aero-' , the excitation energy from bond- 
ing to antibonding molecular orbitals, using the recently 
developed resonant inelastic x-ray scattering (RIXS) 
technique&S RIXS is ideally suited for this study, since it 
provides element-specific and momentum-dependent in- 
formation for electronic excitations^*^ By tuning the inci- 
dent photon energy to the Cu-K absorption edge, one can 
gain information concerning excitations localized around 
the Cu sites without suffering from problems due to the 
La bands. Furthermore, the momentum-resolving capa- 
bility of RIXS provides additional information: the dis- 
persion of such molecular orbital (MO) excitations. In 
this work, we find that Ao-o-* exhibits a strong, systematic 
dependence on the Cu-0 bond length (dcn-o), increas- 
ing as dcu-o is decreased. In addition, for materials with 
a small dcu-o and correspondingly large A^^* , a rela- 
tively large dispersion of the MO excitation is observed. 
We discuss the implication of these observations for un- 
derstanding the electronic structure of the cuprates. 



II. EXPERIMENTS 

The RIXS experiments were carried out at the Ad- 
vanced Photon Source on the undulator beamline 9IDB. 
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TABLE L The copper oxide samples are listed along with the Cu-O bond lengths taken from the references. The top eight 
materials possess perfect square copper oxygen plaquettes, while the bottom three have distorted square plaquettes. Also listed 
are the polarization and the energy of the incident photon in RIXS measurements. 



Label 


Sample 


Crystal" 


dcu-o (A) 


Ref. (dcu-o) 


Polarization 


(eV) 


Ref. (RIXS) 


2122 


Sr2Cu02Cl2 


M 


1.9858 


12 


e ± z 


9001 


This work 


Nd 


Nd2Cu04 




1.9705 


13 


e ± z 


8990 


JA 


Ca 


Ca2Cu02Cl2 




1.9344 




e 1 ~ z ' 


8996 




2342 


Sr2Cu304Cl2 


M 


1.929 


il 


e ± z 


8998 


This work 


LCCO 


Lai.9Cai.iCu206 


B 


1.913 


^8 


e II z 


8999 


This work 


LCO 


La2Cu04 


T 


1.904 


19 


e II z 


8997 


This work 


LSC05 


Lai.95Sro.o5Cu04 


C 


1.898 


i2 


e II z 


8997 


This work 


LSC017 


Lai.83Sro.i7Cu04 


C 


1.885 


J2 


e II z 


8997 


This work'' 


Li 


Li2Cu02 




1.9577 


20 


e ± z 


8997 


21 


CGO 


CuGeOa 




1.9326 


22 


e II x + z 


8990 


23 


112 


SrCu02 


u 


1.910/1.930/1.961 


24 


6 II Z 


8996 


This work 



"The crystals studied in this work were provided by various 
groups, which are denoted here as B:Brookhaven; CiCRIEPI; 
M:M1T; TiToronto; U;Umv. of Tokyo. 

''Since the polarization direction was in the scattering plane in 
this experiment, it changed as momentum transfer was varied. 

<^The RIXS data were taken at T = 15 K. 

'^SrCu02 has three different copper-oxygen bond lengths, which 
are represented as large error bars in Fig. 2. 



Experimental details have been described elsewhere. -W' 
Single crystal samples used in our measurements are 
listed in Table along with several samples studied in 
earlier RIXS experiments. In Tabled c^Cu-o and the ex- 
perimental configuration is listed for each material. All 
measurements were performed at room temperature ex- 
cept for those on the LSC017 sample. In our RIXS ex- 
periments, the scattering plane was vertical and the po- 
larization of the incident x-ray, e, was perpendicular to 
the scattering plane. The polarization direction was kept 
fixed along the direction specified in Table HI where the 
coordinate system reference is the dx^_y2 orbital. That 
is, the copper oxygen plaquette lies in the a;y-plane, while 
the z-direction is perpendicular to the plaquette. We use 
the notation of reduced momentum transfer q through- 
out this paper, with the (tt 0) direction along the Cu-O 
bond direction. 

Before discussing the experimental results in detail, it 
is useful to first review the second-order RIXS process 
to understand the nature of the observed excitation. In 
the ground state, the holes are located in the antibond- 
ing molecular orbital which is a combination of a Cu hole 
state (d^) and an oxygen hgand hole state (d^^L), with 
more weight on the (P state. In the intermediate state of 
this resonance process, a Cu Is electron is excited to the 
Cu 4p band, and the core hole potential alters the balance 
between the and the d}^L states. Then the lowest en- 
ergy state is predominantly \sd}'^LAp, which is lower than 
IsdPAy. These states form the so-called well-screened and 
poorly-screened features, respectively, of the Cu K-edge 
x-ray absorption spectra (XAS). As discussed in detail 
by Hill and coworkers fSiS^ these intermediate states can 
decay into an excited state in which the hole in the anti- 
bonding molecular orbital is filled with an electron, cre- 



ating a hole in the bonding orbital, and an energy loss in 
the outgoing photon. The RIXS process thus creates a 
charge-transfer excitation from bonding to antibonding 
molecular orbitals. 

In our measurements, we have carefully studied the 
incident energy dependence in order to determine the 
resonance energy, for which the MO excitation has the 
maximum intensity. In most cases, the resonance en- 
ergy, which is listed in Tabled corresponds to the higher 
energy peak (poorly-screened feature) in the XAS. On 
the other hand, the intensity of lower energy excita- 
tions near the charge transfer gap (~ 2 eV) is strongly 
enhanced when the incident energy corresponds to the 
well-screened intermediate state, as reported in earlier 
studies." The resonance energy depends on the direc- 
tion of the polarization vector. Detailed results of the 
incident energy dependence study will be published else- 
where. However, one should note that the energy-loss as- 
sociated with the MO excitation does not depend on the 
incident polarization of the photon, although different 4p 
states (e.g., Ap^ or 'Ip-n) are involved in the intermediate 
state, as the polarization is varied with respect to the 
a;y-planei2& 

III. RESULTS 

In Fig.Q] representative RIXS scans are plotted. These 
are energy-loss scans taken at a fixed momentum transfer 
with the incident energy of the x-ray photon fixed at the 
values listed in Tabled The momentum transfer for all 
these scans has been fixed at q=(7r 0), which is the mini- 
mum energy position. The most striking feature in Fig.^ 
is the large shift of the excitation energy from ^ 6 eV for 
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FIG. 1: RIXS spectra taken with the incident energy as spec- 
ified in Table|l]for a fixed reduced wave vector of (vr 0). Each 
spectrum is offset vertically for clarity, and solid lines are 
fits to a Lorentzian lineshape as described in the text. The 
dashed line is a representative scan through the elastic line, 
which shows instrumental energy resolution. 

2122 to - 8 eV for LSC017. To analyze this shift quan- 
titatively, we have fitted the observed excitation spectra 
to a simple Lorentzian lineshape and extracted the peak 
positions, which are plotted in Fig.l^a) as a function of 
Cu-0 bond length, dcu-o- 

A few comments are in order regarding the data anal- 
ysis. First, as is evident from the instrumental resolution 
plotted as a dashed line in Fig.Q the observed excitations 
are not resolution-limited, hence justifying our simple fit- 
ting procedure, which does not convolve the data with 
the instrumental resolution. Second, in several cases, we 
observe more than one type of excitation in these scans. 
For example, the SrCu02 data clearly shows a second 
feature around ~ 8.2 eV, in addition to the main peak 
around 6 eV. In the 2342 case, there also seems to be 
two additional features, one at higher energy (~ 8 eV) 
and the other at lower energy {^-^ 4 eV). However, it is 
difficult to identify these weak features (if present) in the 
data for the other samples, and we have chosen to fit 
all the scans with a single Lorentzian peak with a broad 
width. The peak positions extracted from our analysis 
are, therefore, those of the dominant features. In ad- 
dition, the peak width extracted from the fits might in 
some cases arise from a distribution of several peaks over 
a wide energy range, rather than from a finite inverse life- 
time of a single excitation. Finally, as discussed below, 
we observe dispersion of the MO excitation with momen- 
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Cu-O distance (A) 

FIG. 2: (a) The value of the peak position at (tt 0), Acra* , and 
(b) the amount of dispersion along the Cu-O bond direction 
for each sample is plotted as a function of dcu-o. The solid 
symbols are used for the perfect square plaquettes, while the 
open symbols are for the samples with distorted plaquettes. 
The solid and dashed lines are fits to power law expressions 
Ao-cr* ~ o'cu-o ^^'-^ ^CvL-O' respectively. Note that the mea- 
sured dispersion of the edge-sharing chain compounds Li and 
CGO is not along the Cu-O bond direction, while dispersion 
was not measured for Nd and Caf^*^ 



tum transfer. This can be as large as ~ 0.5 eV in some 
of the compounds studied, as shown in Fig. |31 Thus, the 
peak position of the MO excitation depends not only on 
the sample, but also on q. The scans shown in Fig. ^ all 
are taken at the minimum energy position, q = (tt 0), 
and the peak positions plotted in Fig. IHa) are the values 
measured at this position. 

As shown in Fig. [2Ja) , the energy of the MO excitation 
exhibits a strong dependence on dcu-o- It is notewor- 
thy that the excitation energy exhibits such a system- 
atic dependence on the local structure, and that it is ap- 
parently insensitive to whether the crystal has a planar, 
corner-sharing, or edge-sharing chain strucutre. This is 
consistent with our assignment of these features as MO 
excitations localized within a single Cu-O plaquette. 

The overall trend exhibited in Fig. [2Ja) is not un- 
expected. Intuitively, as the Cu and O atoms move 
closer, the p — d overlap will increase, and the Cu-O 
bonding becomes more covalent with a larger energy 
splitting Aq-o-. . What is surprising is how strong this 
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(271 0) (TtO) (0 0) (JI7I) (27t2jt) 

reduced wavevector q 

FIG. 3: The observed dispersion of the MO excitation. The 
left panel is along the (tt 0) direction, and the right panel 
is data taken along the (tt tt) direction. The solid lines are 
guides to the eye. 



rfcu-o dependence is. We have modeled the dcu-o- 
dependence of Ag-cr* as a power law (Ao-o-* ^ dP) and 
find 77 = —8(2), shown as a solid hne in Fig.[21[a). Note 
that the MO excitation energy is expected to be given 

byS A^^. = y^(2ipp - Ao)^ + IGi^^, where tpp w 0.65 

eV is the hopping matrix element between the oxygen 
p orbitals.^ Since 16^^^ ^ {2tpp — Aq)^, this expres- 
sion leads to Ao-o-* ~ 4^^^, to a first approximation, 
and one then expects similar a d-dependence for Ao-o-* 
and tpd- Our results then implies that the d-dependence 
of tpd is much stronger than that expected from tight- 
binding theory, for which ~ d~^'^ is predicted*?® As plot- 
ted in Fig. El^a), the observed RIXS data clearly devi- 
ates from the Ao-o-- ~ d~^-^ behavior (dashed line). We 
also note that the tpd ~ behavior determined from 
our RIXS measurements is different from the earlier re- 
port of tpd ^ d~'^ by Cooper and coworkers, which was 
estimated indirectly from a three-band Hubbard model 

2Q 

expressionisa 

In Fig. 131 the dispersion of the MO excitation is plotted 
as a function of q for selected samples. These data sug- 
gest that the picture of a completely localized MO excita- 
tion is an oversimplification - such a model would predict 
no dispersion of these features. The data of Fig. [2Ib), 
which show significant dispersion for dcu-o < 1-93 A, 
suggest then that this localized picture breaks down for 
small bond distances. For the LSC017 sample, the dis- 
persion bandwidth is about 0.5 eV, with the minimum 
excitation energy occuring at the zone boundary, imply- 



ing an indirect gap. Note that this is completely different 
from the direct nature of the lower energy charge-transfer 
gap as reported in Ref.^J As A^a* decreases (e.g., LCO 
and LCCO), the dispersion of the MO excitation becomes 
weaker. For the samples with an even smaller , the 
dispersion becomes flat, as shown for the 2342 sample 
(Figs. 2 and 3). Limited momentum dependence mea- 
surements for the 2122 sample (not shown) also show a 
dispersionless behavior. The size of the observed disper- 
sion along the Cu-0 bond direction is plotted against 
dcu-o in Fig.Etb). 

As shown in Fig. El the size of the dispersion and Ao-o-* 
both increase as dcu-o decreases. This suggests that the 
bandwidth of the dispersion is also controlled by the hop- 
ping parameter tpd- Such behavior would, of course, be 
expected from the increased overlap of the wavefunctions, 
since charge carriers then become less localized. How- 
ever, the simplest picture of delocalized electrons fails 
to describe the observed dispersion. For example, the 
bandwidths of the Cu-0 bonding and antibonding bands 
in LCO are very large 3 eV, and interband transitions 
between these two bands would have a direct gap of ~ 4 
eV^ 



IV. DISCUSSION AND SUMMARY 

One of the most surprising results in this study is 
that the l^aa' — c^Cu-o scaling seems to apply to dif- 
ferent structures. In contrast, previous studies of the 
bond-length scaling of various quantities, such as charge- 
transfer gap (Apr), or superexchange interaction (J), 
have been limited to compounds with corner-sharing 
structures.^ It is well known that such quantities as su- 
perexchange coupling depend not only on the p — d hy- 
bridization, but also crucially on the angle between the 
two Cu-0 bondsi^ One can argue that A„„* is a better 
measure of Cu-0 hybridization than J or Act, since it 
is only dependent on dcu-o Sind not on the presence of 
neighboring atoms. 

We have noted that the tpd ~ dependence in- 
ferred from our study deviates significantly from the 
tight-binding picture. It also appears to give rise to 
discrepancies with other experiments. For example, for 
materials with the corner-sharing structure, the three- 
band Hubbard model gives J ^ f^^/U* H^^jj,, where U* 
is an effective onsite Colomb interaction which is as- 
sumed to be constant. If we use the experimentally 
determined^S Act ^ d^^ and our tpd ~ d~^ result, 
we obtain J ^ d~^^, which clearly disagrees with the 
much weaker d-dependence of J observed in various ex- 
periments, including two magnon Raman scattering^ 
One might expect that this discrepancy could be re- 
solved by considering the fact that, due to strong elec- 
tron correlations, the simple tight-binding picture of co- 
valent bonding has to be modified. In fact, Mizuno and 
coworkers'^'^ considered two contributions to tpd- That 
is, tpd — -I- ip^, where is the contribution from the 
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atomic potential which depends only on dcu-Oj while 
is the contribution from the Madelung potential, which 
depends on the detailed arrangement of the neighbor- 
ing ions. However, the calculated contribution from the 
Madelung potential is of order of ^ 0.1 cV or smaller, '^■^ 
so that this alone is not enough to explain the ^ d^^ 
dependence. 

These results may be suggesting that one has to aban- 
don the simple relationship of Ao-o-* ~ 4ipd. Certainly, 
as discussed above, the picture of a completely local- 
ized MO excitation is apparently an oversimplification, 
since it breaks down as c?cu-o becomes shorter - as ev- 
idenced by the sizable dispersion observed in LSC017. 
Thus, if a more realistic expression for Ao-o-* is used, 
tpd ~ d~^'^ scaling law might be recovered. For exam- 
ple, a recent first-principles calculation has emphasized 
the role of apical oxygens in the systematics of high tem- 
perature superconductivity'^^. Indeed, the scaling plot 
in Fig. 121 also exhibits some systematic dependence on 
the number of apical oxygenS)^ and it may be interest- 
ing to further investigate the role played by apical oxy- 
gens. Certainly, a systematic ab initio calculation of the 
dcu-o-dependence of MO excitation energy in large clus- 
ters would be highly desirable and may help to clarify the 
relationship between Ao-o-* , tpd, and dcu-o- 

To summarize, we have studied a charge-transfer exci- 



tation in various cuprate compounds using resonant in- 
elastic x-ray scattering technique. We assign this exci- 
tation to a mostly localized molecular orbital excitation, 
that is, an excitation from a bonding to an antibonding 
molecular orbital. We have found that this molecular 
orbital excitation energy, which is a measure of the hop- 
ping matrix element tpd, exhibits a systematic Cu-0 bond 
length dependence, which is much stronger than that ex- 
pected from tight-binding theory. We have also observed 
a sizable dispersion of this excitation in some materials, 
suggesting that this molecular orbital excitation becomes 
less localized as the p — d hybridization becomes large. 
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